Neutrino spectra from stellar electron capture 
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Using the recent shell model evaluation of stellar weak interaction rates we have calculated the 
neutrino spectra arising from electron capture on p/-shell nuclei under presupernova conditions. 
We present a simple parametrization of the spectra which allows for an easy implementation into 
collapse simulations. We discuss that the explicit consideration of thermal ensembles in the parent 
nucleus broadens the neutrino spectra and results in larger average neutrino energies. The capture 
rates and neutrino spectra can be easily modified to account for phase space blocking by neutrinos 
which becomes increasingly important during the final stellar collapse. 
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PACS numbers: 26.50.+X, 23.40.-s, 21.60.Cs 

Weak interactions play a central role in the final evo- 
lution of massive stars These processes, mainly elec- 
tron capture and beta decays, create neutrinos which for 
all densities until the collapse becomes truly hydrody- 
namic (i.e. p ~ 10 11 g cm~ 3 ) escape the star carrying 
away energy and reducing the entropy. Thus, in simu- 
lations which follow the star's evolution until the iron 
core reaches central densities of order a few 10 9 g cm -3 
it is sufficient to know the average neutrino energies of 
the various weak processes to determine the energy loss 
rate along the stellar trajectory, e.g. (|]. Such simula- 
tions define the 'presupernova models', e.g. which are 
then being used as input to detailed studies of the col- 
lapse and explosion mechanism, e.g. j| . In recent years it 
has become apparent that the various neutrino reactions 
play an essential role in these simulations and, as neu- 
trino reactions depend sensitively on energy, these studies 
require a detailed bookkeeping of the neutrino spectra. 
This goal is achieved within hydrodynamical models with 
explicit neutrino Boltzmann transport ^-R]. The neces- 
sary input into such simulations are the weak rates and 
the associated neutrino spectra. 

So far, electron capture on nuclei has only been treated 
quite schematically in collapse simulations with neutrino 
Boltzmann transport, e.g. The ensemble of nuclei 

in the stellar composition, given by nuclear statistical 
equilibrium, is represented by an 'average nucleus', whose 
capture rate then is derived on the basis of the indepen- 
dent particle model (even reduced to a model which only 
considers / 7 / 2 and / 5 / 2 orbitals jjjj). The correspond- 
ing neutrino spectrum is then approximated in the spirit 
of this model, assuming a nucleus-independent energy 
splitting of 3 MeV between the / 7 / 2 and / 5 / 2 orbitals. 
We note, that the standard for the stellar weak interac- 
tion rates in presupernova evolutions, however, was set 
by the work of Fuller, Fowler and Newman (FFN) ||| 
who estimated the rates for the nuclei with mass num- 
bers A = 21 — 60 on the basis of the independent par- 
ticle model and experimental data, whenever available. 
Presupernova evolution studies then considered the FFN 
rate tables for a proper nuclear composition. FFN list 
the average neutrino energies for the various weak reac- 



tions, but to our knowledge neutrino spectra have not 
been derived on the basis of the FFN compilation. As 
mentioned above, such spectra are also not required for 
presupernova studies as the neutrinos leave the star un- 
hindered. 

Due to progress in nuclear many-body modeling and 
in computer hardware and guided by experimental data 
it has recently been possible to treat the nuclear struc- 
ture problem involved in the calculation of stellar weak 
interaction rates in a reliable way 0,0. The calcu- 
lations have been performed on the basis of large-scale 
shell model studies which reproduce all experimentally 
available relevant data quite accurately ]l2| ]. These shell 
model rates show some marked differences to the FFN 
estimates, leading to significant ch ang es in the presuper- 
nova evolution of massive stars |@,[L3|. It appears there- 
fore reasonable that this compilation [jllj should also be 
used in collapse and explosion studies which build on 
the presupernova models. To make such use possible we 
will here study the neutrino spectra corresponding to the 
shell model rates and show a way how these spectra can 
be easily and consistently implemented in collapse codes. 
It is important to note that under the conditions of the 
presupernova models and in the subsequent stellar evo- 
lution beta decay is strongly blocked by the appreciable 
electron chemical potential and the total electron capture 
rates are orders of magnitude larger |p^,p|JT^|. Thus it is 
quite sufficient to focus on the neutrino spectra arising 
from electron captures for post-presupernova simulations 
where detailed neutrino transport is important. 

Under the stellar conditions we are concerned with 
electron capture is dominated by Gamow- Teller (GT) 
transitions [ 14j . The appropriate formalism has been de- 
rived in H|) | : 

A = ln2 ^ 2 ^ (2 J, + l)e-^ fcT ) 
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where the sums in i and j run over states in the parent 
and daughter nuclei, respectively. For the constant K 
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we used K — 61 46 ±6 s [|L5| and gv,9A are the vec- 
tor and axialvector coupling constants. G(Z, A, T) = 
^2 li exp(—Ei/(kT)) is the partition function of the parent 
nucleus. The sum in the GT matrix element runs over 
all nucleons. In the phase space integral w is the total 
energy of the electron in units of m e c 2 , and p = ^/w 2 — 1 
is its momentum in units of m e c. Finally the Q- value for 
a transition between two nuclear states i, f is defined in 
units of m e c 2 as 
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where M p , Md are the nuclear masses of the parent and 
daughter nucleus, respectively, while Ei,Ej are the exci- 
tation energies of the initial and final states. The lower 
integral limit in (1) is Wi = 1 if Qy > —1 or wi — \Qij\ 
if Qij < — 1. For the stellar conditions we are interested 
in, the electrons are well described by Fermi-Dirac (FD) 
distributions, with temperature T and chemical potential 



Se(E e ) 



1 



exp 



kT 
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with E e = wm e c 2 . The remaining factor appearing in 
the phase space integrals is the Fermi function, F(Z,w), 
that corrects the phase space integral for the Coulomb 
distortion of the electron wave function near the nucleus. 

Applying the formalism described above and determin- 
ing the nuclear matrix elements within large-scale shell 
model calculations, stellar electron capture rates have re- 
cently been calculated for p/-shell nuclei. These nuclei 
dominate the weak processes in the presupernova evolu- 
tion of massive stars and details of the rate evaluations 
can be found in . Here we will use the same approach 
to study the neutrino spectra emerging in the electron 
capture reactions. Energy conservation requires that the 
neutrino emitted after capture of an electron with energy 
w on an initial state i leading to a final state j is: 



E v = m e c 2 (w + Q^). 



(4) 



The neutrino spectra for a specific nuclear transition 
i — > j is given by the respective partial rate per energy 
interval. The total spectrum is then the sum over all 
possible transitions. The normalized neutrino spectrum 
n(E v ) is obtained by dividing by the total electron cap- 
ture rate. 

In the derivation above we have explicitly assumed that 
the neutrinos produced by the electron capture process 
leave the star; i.e. there is no blocking of the phase space 
due to the presence of neutrinos in the stellar environ- 
ment. The spectra, which we will discuss below, are all 
derived based on this assumption. In the collapse phase 
following the presupernova evolution neutrinos are get- 
ting increasingly trapped in the core. This will require 
the inclusion of a neutrino blocking factor (1 — Si/(E U )) in 



the phase space integral of the electron capture formula, 
Eq. (1). The neutrino distribution S U (E W ) depends on 
position and time and can be calculated within the Boltz- 
mann transport formalism. Once S V (E V ) is known, the 
present neutrino spectra and the capture rates of |Tc| ] 
can be easily corrected for neutrino blocking. For the 
neutrino spectra this is achieved by folding the uncor- 
rected spectra n(E u ), as presented in this paper, with 
the blocking factor (1 — S V (E V )). The corrected capture 
rates are obtained by multiplying the tabulated rates plj] 
with / n(E)(l - S v {E))dE. 

For the following discussion it is useful to realize that 
the neutrino spectra depend basically on three quantities 



• the electron chemical potential 

• the 'effective' Q-value, Q c ff = 

• the GT strength distribution 



M p -M d + Ei 



Generally one expects that neutrino energies are in- 
creased for larger electron chemical potentials, favor- 
able mass differences between daughter and parent nuclei 
(M p — Md), and from thermally excited states. Further- 
more, large neutrino energies are favored if the strong 
GT transition resides at low excitation energies in the 
daughter nucleus. Due to nuclear pairing structure ar- 
guments |l(J this is the case for odd-odd daughter nuclei 
(capture on even-even parents) , while the bulk of the GT 
strength is somewhat higher (~ 2 — 3 MeV) for odd-A 
nuclei and is shifted by additional 2-3 MeV in even-even 
daughter (capture on odd-odd parents) p2[ . We note 
that these differences in the GT distributions result in 
the fact that low-lying rather weak GT transitions con- 
tribute relatively more to the electron capture rates on 
odd-odd nuclei and odd-A nuclei, while for even-even par- 
ents the GT bulk often resides at such low energies (e.g. 
Ni isotopes) that a distiction between low-lying strength 
and GT bulk is not very meaningful. 

To illustrate the discussion we plot in Fig. 1 the 
normalized neutrino spectra for electron capture on the 
A=56 isobars. The calculation has been performed for 
typical conditions during silicon shell burning of a 15 M Q 
star [| (T = 4 x 10 9 K, p = 3 x 10 8 g cm" 3 , Y e = 0.45). 
The resulting electron chemical potential then is p e = 2.5 
MeV. 

We note that 56 V, 56 Cr, 56 Mn and 56 Fe have Q Q = 
M p — Md < 0, thus making electron capture in the labo- 
ratory impossible. The other two nuclei 56 Co (Qq — 4.06 
MeV) and 56 Ni (Q = 1.62 MeV) decay dominantly by 
electron capture. With the exception of 56 Co the neu- 
trino spectra for the other 5 nuclei are very similar: 
they are peaked around rather small neutrino energies 
E„ « 1-2 MeV with a width of 1.4-1.8 MeV. The reason 
for this quite similar structure is twofold. For the nuclei 
with negative Qo-values, electron capture is hindered and 
requires electrons from the exponentially decreasing tail 
of the FD distribution. Obviously to achieve an apprecia- 
ble rate it is advantegeous to keep the neutrino energies 
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FIG. 1. Normalized neutrino spectra for stellar electron 
capture on selected A — 56 isobars. The spectra have been 
calculated for stellar conditions which are typical for silicon 
shell burning in a 15 Mq star P|. 



small. Thus capture to low-lying states occurs with elec- 
trons with lower energies than capture to the bulk of the 
GT strength, but both are accompanied by low-energy 
neutrinos. For 56 Ni the Qcr value allows capture of elec- 
trons with all energies. However, the GT distribution 
in the daughter 56 Co is well concentrated at low excita- 
tion energies, resulting again in a rather narrow neutrino 
spectrum. 
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FIG. 2. Partial contributions of individual states in the 
parent nucleus to the neutrino spectrum for stellar electron 
capture on 56 Co. The calculation has been performed for the 
same conditions as in Fig. 1. The spectra are multiplied by 
their relative weight to the total capture rate. 

The 56 Co spectrum is quite different, showing a double- 
bump structure. The reason is explained in Fig. 2 
which shows the partial neutrino spectra contributed by 
selected states in the parent nucleus. The spectra for 
both, the ground state (J = 4) and the first excited 
state (J = 3) show the single-peak structure, as gen- 
erally observed for the other nuclei. We also find the 
spectra for the first excited state somewhat wider than 



for the ground state. This reflects the strong angular mo- 
mentum mismatch for the ground state which does not 
connect to states in 56 Fe below 2 MeV and has only very 
weak transitions to states below 3 MeV. This is differ- 
ent for the excited state which can reach the lower-lying 
J = 2 states in 56 Fe and in that way produces neutri- 
nos with larger energies. Interestingly the excited J = 1 
state at 1.7 MeV excitation energy in 56 Co produces a 
neutrino spectrum with a peak energy around E v = 8 
MeV. We note that this state has a rather strong GT 
matrix element to the 56 Fe ground state and hence, con- 
sidering its favorable effective Q value of Q e g ~ 5.8 MeV, 
electron capture on this state generates neutrinos with 
rather large energies. While the excitation energy in the 
parent increases the effective Q-value, and thus the av- 
erage neutrino energy, it reduces the contribution to the 
rate due to the Boltzmann weight. However, the excited 
1 + state yields the clue to the higher-energy neutrino 
peak. We notice that there are many more states in the 
excitation energy range ~ 2-4 MeV which are connected 
to the low-lying states in the daughter nucleus 56 Fe by 
strong GT transitions. FFN have coined the term 'back- 
resonances' for these states [||9| as they are part of the 
bulk of the GT strength built on the low-lying states in 
the inverse direction (Ref. |l^] explains how these states 
are considered in the rate evaluation.) Electron capture 
on these backresonances occurs with favorable Q c ff- value 
and hence allows the emission of neutrinos with rather 
high energies. Despite the Boltzmann suppression the 
gain in phase space combined with the large matrix ele- 
ments ensure that these states combined contribute no- 
ticeably to the total rate and produce the second peak 
in the total neutrino spectrum. Noting that the inte- 
gral over the spectrum reflects the relative contribution 
to the total rate, we remark in passing that Fig. 2 also 
implies that at the chosen conditions electron capture on 
56 Co is dominated by the one on the ground state. Al- 
though similar in nuclear structure, the electron capture 
on the odd-odd nuclei 56 Mn and 56 V does not produce 
a double-bump structure due to the negative Qq value 
which favors emission of low-energy neutrinos. 

At the conditions of silicon shell burning, depicted in 
Fig. 2, the electron chemical potential is yet not large 
enough (fi e — 2.5 MeV) to allow significant capture 
on the odd-odd nucleus 56 Co from low-lying states to 
the bulk of the GT + distribution in the daughter nu- 
cleus which resides at around ~ 7-9 MeV in 56 Fe [ jl6| . 
However, in the subsequent stellar evolution /x e increases 
rather fast. Thus, capture to the bulk, e.g. with signif- 
icantly larger GT strength, becomes easier for the low- 
lying states. This increases their relative weight com- 
pared with the one of the backresonances and is also not 
compensated by the relative gain in the Boltzmann factor 
of the latter. This behavior is demonstrated in Fig. 3, 
again for 56 Co. The temperatures have been chosen ac- 
cordingly, using the stellar trajectories as given by Heger 
et al. Q. We also observe that, once capture to the 
bulk of the GT distribution dominates the capture rate, 
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FIG. 3. Normalized neutrino spectra for stellar electron 
capture on 56 Co at several different phases of the final evo- 
lution of a 15 Mq star. The stellar parameters have been 
taken from Table 1 of The chemical potentials are given 
in MeV, while Tg defines the temperature in 10 9 K. 

the spectrum becomes single-peaked (approximately for 
fi e — 4.2). A further increase in p e then simply increases 
the average neutrino energy (e.g. for p e = 8.1). 

We summarize that the relative height of the two peaks 
reflects the competition between electron chemical po- 
tential (its decrease reduces the capture from low-lying 
states to the bulk of the GT distribution) and tem- 
perature (its decrease reduces the Boltzmann weight of 
the backresonances) . If one follows the stellar evolution 
backwards in time, i.e. to smaller temperatures, den- 
sities and electron chemical potential, the spectrum ul- 
timatively becomes single-peaked as it is dominated by 
the temperature-favored ground state contribution. The 
spectrum also develops 'discontinuities' (see e.g. for p, e = 
0.9) which reflect the fact that the spectrum represents a 
sum over several initial and final states which all have a 
definite minimal neutrino energy £™ m = m e c 2 (Qij + 1). 
At -E™ 111 , the electron momentum vanishes (p — 0), but 
pF(z,u) is finite resulting in a finite value for n(E 1 J Lln ). 
For larger p, e values these discontinuities are smeared out 
as the individual spectra noticeably overlap. 

As examplified for 56 Fe in Fig. 4 the situation is quite 
different if one studies the neutrino spectrum emerging 
from capture on an even-even nucleus as function of elec- 
tron chemical potential. To understand the reason we 
note two facts. First, as the GT distribution in the 
daughter is quite concentrated at low excitation energies 
a double-peak structure related to the distinct capture to 
low-lying states and the GT bulk does not emerge. Sec- 
ond, in even-even nuclei the backresonances are at higher 
excitation energies than in odd-odd nuclei |l(| . The con- 
sequences for the neutrino spectrum are obvious. If com- 
pared to 56 Co, the relative contribution of the low-lying 
states to the electron capture on 56 Fe is significantly en- 
hanced with respect to the backresonances and no high- 
energy neutrino peak emerges. As an increase of the 
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FIG. 4. Normalized neutrino spectra for stellar electron 
capture on 56 Fe at several different phases of the final evolu- 
tion of a 15 M0 star. The stellar parameters have been taken 
from Table 1 of Q . The last panel corresponds to typical con- 
ditions during the collapse phase (T = 10 10 K, p = 3 x 10 10 
g/cm 3 , and Y e = 0.42). The chemical potentials are in MeV, 
while Tg defines the temperature in 10 9 K. 

electron chemical potential makes the capture energeti- 
cally easier, the neutrino peak energy moves to higher 
energies. We mention that the last panel of the figure 
already corresponds to a phase of the contraction after 
the presupernova model which we have approximated by 
T = 10 10 K, p = 3 x 10 10 g/cm 3 and Y e = 0.42. It 
demonstrates that the simple structure of the spectrum 
remains also during that stellar evolution stage which re- 
quires detailed neutrino transport. 

As stated above detailed neutrino transport becomes 
important in the final evolution of massive stars, follow- 
ing the presupernova models. Due to electron captures 
the matter in the final presupernova models is neutron- 
rich and the nuclei present have negative Qo-values. On 
the other hand, the electron chemical potential has grown 
strongly enough until this point thus allowing capture to 
the bulk of the GT strength. Fig. 5 shows the neu- 
trino spectra for the 6 nuclei which dominate electron 
capture in the presupernova models of a 15 M star M. 
Due to Heger et al. the core density and temperature 
are p = 9.1 X 10 9 g/cm 3 and T = 7.2 x 10 9 K, while 
the Y e value is 0.432 j|. The chemical potential then 
is p e — 8.1 MeV. We note that all neutrino spectra are 
single-peaked. The average neutrino energy released by 
nuclei is about 3 MeV, while it is 6.25 MeV for capture on 
free protons which under these presupernova conditions 
become abundant enough to significantly contribute to 
electron capture. 

As neutrino cross sections scale with E%, high-energy 
neutrinos are more easily trapped. Capture on free pro- 
tons has a more favorable Qo value than capture on neu- 
tronrich nuclei present in the presupernova matter com- 
position. As a consequence the neutrino average energy is 
higher for capture on free protons. Nevertheless captures 
on nuclei still produces the larger amount of high-energy 
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FIG. 5. Normalized neutrino spectra for stellar electron 
capture on the six most important 'electron-capturing nu- 
clei' in the presupernova model of a 15 M© star, as iden- 
tified in |). The stellar parameters are T = 7.2 x 10 9 K, 
p = 9.1 X 10 9 g/cm 3 , and Y e = 0.43. The solid lines repre- 
sent the spectra derived from the shell model electron capture 
rates. The dashed line shows the fit to the spectra, using the 
parametrization of Eq. (6) and adjusting the parameter q to 
the average neutrino energy of the shell model spectrum. The 
dashed-dotted spectrum corresponds to the parametrization 
recommended in [Hi. 



neutrinos in lighter stars (M < 20 M ) as they dominate 
the rate in the presupernova models. This is different for 
heavier stars where the capture on free protons in the 
presupernova model corresponds already to 30-50 % . 

For presupernova conditions the shell model evaluation 
predicts slightly larger average neutrino energies than the 
FFN rates A possible explanation for this difference 
is given by the fact that the shell model rates explicitly 
consider capture from thermally excited states. The Q- 
value for 59 Fe is Qq — —5.696 MeV. However, the excited 
states have more favorable Q c s values than the ground 
state thus making electron capture more easy and in sev- 
eral cases supporting larger neutrino energies. Fig. 6 
shows the neutrino spectrum calculated for the 6 lowest 
states in 59 Fe. (Here the distributions are multiplied by 
their relative weight in the total rate). One finds that 
indeed neutrino spectra from excited states often have a 
wider tail. This results from capture to low-lying states 
in the daughter. In contrast to the low- lying transitions 
the GT bulk approximately obeys Brink's hypothesis, 
e.g. |Tlj| , This states that the GT strength of excited 
states is the same as for the ground state, only shifted by 
the excitation energy of the parent state. In particular, 
Brink's hypothesis implies that the relevant energy dif- 
ference Qeff — Ej is the same for capture to the GT bulk 
for all parent states. Another observation has already 
been mentioned above. Although the total capture rate 
is still dominated by the transition to low-lying states, 
the main source for high-energy neutrinos, however, are 
the 'backresonances' pB. 

If the shell model neutrino spectra are to be used 
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FIG. 6. Partial contributions of individual states in the 
parent nucleus to the Neutrino spectrum for stellar electron 
capture on 59 Fe. The calculation has been performed for the 
presupernova conditions of Fig. 5. The spectra are multiplied 
by their relative weight to the total capture rate. The left 
panel shows the neutrino spectra calculated for the ground 
state (solid line) and the 5 lowest excited states in 59 Fe. The 
right panel compares the total neutrino spectrum (solid line) 
with the one obtained only from the backresonances (dotted 
line, for a definition see text and Jl(|). The dashed line shows 
the contributions from the low-lying individual states which 
have been considered to calculate the electron capture rate 

in collapse simulations, they have to be represented by 
parametrizations which are accurate, fast and can be eas- 
ily implemented. Our proposal for the parametrization 
is based on the following approximations. Suppose that 
i) the electron capture on a state in the parent nucleus, 
described by GT transitions, leads to a single state in 
the daughter nucleus at energy Ej and that ii) Brink's 
hypothesis is valid, i.e. this state is at Ej + Ei if the 
capture is on an excited state in the parent at excitation 
energy Ei. If we equal electron energy and momentum, 
which is a valid approximation for the conditions we are 
interested in, the neutrino spectrum has the form Q 



n{E v ) = El{E v - qf 



N 



l + exp{(£„ -q~ii e )/kT} 



(5) 



with q = Qq — Ej and a constant N which normalizes 
the neutrino spectrum to unity. This form is obviously 
valid for capture on free protons where the parameter q 
is the reaction Qq value, q — —1.29 MeV. For finite nu- 
clei, q should be considered a fit parameter. It can be 
adjusted to the average neutrino energy which is listed 
in the shell model rate tabulations for a grid of tem- 
perature/density points and can be easily interpolated 
in-between. We have tested this proposal and gener- 
ally find that the parametrization approximates the shell 
model spectra rather well, as can be seen by the dashed 
curves in Fig. 5. Of course, our parametrization fails if 
the spectrum is double-peaked as observed for capture 
on odd-odd nuclei under special conditions (e.g. see Fig. 
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1.). However, we do not expect that these conditions oc- 
cur for p/-shell nuclei during the collapse phase where 
the electron chemical potential is high enough compared 
to the reaction Qo value to allow appreciable electron 
capture to the GT bulk. Finally we note that for 52 V the 
shell model spectrum is wider than the parametrization. 
This is caused by the fact that the spectrum for the in- 
dividual states shows noticeable differences and does not 
strictly follow the Brink hypothesis. For example, the 
excited state of 52 V at 22 keV has J = 5 and thus there 
are no low-lying states in the daughter 52 Ti which can 
be reached by GT transitions. This is different for the 
excited J = 1 state at 141 keV which connects strongly 
to the 52 Ti ground state. 

Bruenn suggested a similar parametrization for the 
neutrino spectra emerging from electron capture on pf- 
shell nuclei, however simply setting q = Qo — 3 MeV ||. 
The resulting spectra are compared to the shell model 
spectra in Fig. 5. Despite the simple guess for the pa- 
rameter, the agreement is quite acceptable. 

In summary, the knowledge of the neutrino energy 
spectra at every point and time in the core is quite rel- 
evant for simulations of the final collapse and explosion 
phase of a massive star. In the collapse phase, neutri- 
nos are mainly produced by electron capture on nuclei 
and protons and their emerging energy spectra are an 
important ingredient in the simulations. In this paper 
we have presented neutrino spectra for stellar electron 
capture during the final presupernova evolution stage of 
massive stars. The spectra have been consistently derived 
in the framework of the recently evaluated capture rates 
which have been calculated on the basis of state-of-the- 
art large-scale shell model studies. Furthermore we have 
calculated the spectra for stellar conditions which have 
been obtained in presupernova evolution of massive stars, 
using the same shell-model weak interaction rates. The 
calculated presupernova neutrino spectra show a rather 
simple structure which is easily parametrizable. This 
parametrization is easily implementable into the simu- 
lation codes and allows for a derivation of the neutrino 
spectra consistent with the shell-model weak-interaction 
rates. 
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